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DEFINITION OF SYMBOLS 


Symbol 

Definition 

Units 

a 

surface area to unit volume ratio of screen 
wire (defined in Ref. 2) 

cm -1 

A 

area 

cm 2 

A' 

actual area open to flow, Ae 

cm 2 

A l» a 2 

empirically determined constants (defined in 
Ref. 4) 

dimensionless 

Ac 

constant, A^L/D a 2 g c (defined in Ref. 4) 

cm -1 

B 

screen thickness (defined in Ref. 2) 

cm 

B c 

constant, A 2 L/D a g^ (defined in Ref. 4) 

dimensionless 

D 

depth of channel 

cm 


characteristic pore size 

cm 

D h 

hydraulic radius (defined in Ref. 9) 

cm 

Eu 

Euler number, 2Ap/pV 2 

dimensionless 

f 

friction factor Apg c £ 2 D/TBpV 2 (defined in 
Ref. 2) 

dimensionless 

F 

parameter, fH/4D^ (defined in Ref. 9) 

dimensionless 

g 

acceleration due to gravity 

cm/ sec 2 

Sc 

lbmft 

conversion constant, e.g., 32.174 , 

Ibfsec 2 

depends on system 
of units used 

H 

channel length (defined in Ref. 9) 

cm 

K 

constant, 2LAi/D a (defined in Ref. 8) 

dimensionless 

Ko 

screen parameter, Ap/V 

gm/ cm 2 sec 

L 

channel width (definition in Ref. 9) 

cm 


ii 


Symbol 

°P 

Ap 

Q 

<T 

Re 

Re a 

S 

T 

u 

V 

V e 

x,y 

z 


Definition 

fraction of area open to flow (defined in 
Ref. 1) 

screen pressure drop 

mass flow rate (defined in Ref. 9) 

volumetric flow rate 

pV 

Reynolds number, " ~ 2 (defined in Ref. 2) 

|J3. U 

Reynolds number, D a Vp/y 

solidity, 1 - Op (defined in Ref.l) 

tortuosity factor (defined in Ref. 2) 

fluid velocity in x-direction 

screen approach velocity — 

A 

entering screen velocity, iL. 

Ae 

Cartesian-coordinate distances 
length, X/H (defined in Ref. 9) 


a 

e 

e 

n 

y 

v 

P 

0 


viscous coefficient (defined in Ref. 2) 

inertial coefficient (defined in Ref. 2) 

screen void fraction (defined in Ref. 2) 

constant 

fluid viscosity 

fluid kinematic viscosity 

fluid density 

parameter, 2QH/D Z LK 0 (defined in Ref. 9) 


Units 

dimensionless 

gm/ cm. sec 2 

gm/sec 

cm 3 / sec 

dimensionless 

dimensionless 

dimensionless 

dimensionless 

cra/sec 

cm/ sec 

cm/sec 

cm 

dimensionless 

dimensionless 

dimensionless 

dimensionless 

dimensionless 

centipose 

cm 2 /sec 

gm/cm 3 

dimensionless 


iii 


Symbol 


Def initlon 


* 

GDCD 

MDAC 


denotes average value when placed above 
symbol, e.g. V is average velocity 

denotes dimensionless quantity, e.g. V* is the dimensionless 
velocity. 

Convair Division of General Dynamics 
McDonnell Douglas Astronautics Company 


NRSD 


North American Rockwell Space Division 


CHAPTER I 


INTRODUCTION 


The purpose of the present investigation is to gather experimental 
pressure drop and velocity data for woven screens. 

Previous investigators have presented correlations for predicting 
pressure drop across woven screens. All correlations presented to date 
have been developed from "pipe flow" (i.e., circular) configurations and 
data taken in those experiments. The present experimental investigation 
uses a rectangular channel configuration. 

The report is divided into two major sections and an appendix A and B. 

Section I presents experimental data taken for three dutch twill 
screens (50 x 250, 200 x 600, and 325x 2300) and two square weave screens 
(200 x 200, and 325x325) using tap water as the test liquid. Pressure drop 
measurements were made at four locations in a rectangular channel 8.89 cm 
long. In every case, no variation (for a given screen) in pressure drop 
at the four locations was measured. The data is presented as Ap versus 
V„ where V is the average entering velocity and is calculated by dividing 
the volumetric flow rate by the screen area open to flow. Existing data 
have been based upon an average approach velocity (volumetric flow rate 
divided by total screen area) . 

Section II presents experimental data for a 50 x 250 dutch twill 
screen using water in a modified experimental apparatus. The channel length 
is extended to a length of 29.16 cm to study the effect of pressure drop 
variation as a function of Z (dimensionless length along the screen) , 
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Channel depth is made variable in order to study its effect upon pressure 
drop. Basic data is presented as Ap versus Z. 

The equations of continuity and momentum for the present experi- 
mental model are presented in Appendix A. Also included in Appendix B, 
is a computer program listing of an extension of a McDonnell Douglas 


theoretical model and data from that computer program. 



SECTION 1 


PREDICTION OF PRESSURE DROP ACROSS WOVEN SCREENS 


1-1 Previous Investigations 

Several investigators have attempted to obtain empirical or analytical 
relations between the screen pressure drop and velocity. The correlations 
are based upon an average approach velocity normal to the screen (i.e., the 
volumetric flow rate divided by the total screen area) and screen parameters. 
The empirical equations below are predominant in the literature. 


£ - £ + l? 


(l-D 


where f is the friction factor defined by a P §c £ p 


T B p V 


PV 


Re is the Reynolds number defined by — 5 - 

pa z D 


a, 3 are constants determined experimentally 


eu = 

PV 2 


(1-2) 


where Eu is the Euler number and is defined in the following manner 
for square weave and dutch twill screens 


Eu = 


IWI 

to 

II 

7 

03 

1 

1 

0. 

0 

1 


Li-sj 


for square weave screens. 


(1-3) 


Equation (1-3) is the Hoerner equation (Ref. 1) where Op is the 
fraction of area open to flow and S is the solidity. 
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Eu 


8L 

dT 


- — + k 
Re a 2 


for dutch twill screens 


(1-4) 


where 


L is the thickness of the screen 

D a is the characteristic pore size 

A]_, A 2 are constants determined experimentally 

Re a is the Reynolds number based on D a and defined as D a Vp/y. 

Equation (1-1) was developed by Armour and Cannon (Ref. 2) who modeled 
screens as a collection of submerged objects with surface area to unit volume 
ratio "a" for laminar flow and as a bundle of tubes of diameter D for 
turbulent flow. Pressure drop data for the flow of gaseous nitrogen and 
helium through plain square, full twill, fourdrinier, plain dutch, and dutch 
twill screens were used to derive coefficients of a = 8.61 and 3 = 0.52. An 
illustration of the types of screens is shown in Eigure (1-1) . 

Other investigators have arranged their data in the form of the 
Armour-Cannon correlation. McDonnell Douglas Astronautics Corporation 
(MDAC) , using GN 2 and He as test fluids, presented data on three dutch 
twill screens as shown in Figure (1-2) (Ref. 3)*. Their data points for the 
friction factor, f, were lower than the Armour Cannon correlation (generally 
by a factor of 2 for each Reynolds number; however, the correlation was 
successful in aligning the data points for the three dutch twill screens. 


* Ref. 3, Fig. 33, p. 76. 
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General Dynamics Convair Division (GDCD) tested the following six 
dutch twill screens in 1969 (Ref. 4)' 

• 80 x 700 

• 165 x 800 

• 150 x 700 

• 30 x 250 

• 50 x 250 

• 20 x 250 

using GH^ and GN^ as test fluids. They compared the results of those tests 
with previous test data taken in 1968 on three other dutch twill screens 
and one square weave screen (Ref. 5): 

• 200 x 1400 

• 165 x 1400 (Dutch Screen) 

• 200 x 600 

• 20 x 20 (Square Screen) 

using GHe, LH2, GI^, LID,, G^, tap water, and distilled water as test fluids. 
Their results of friction factor, f, as shown in Figure (1-3)* were also 
bilow the Armour-Cannon correlation, and they found that the data could best 
be represented by the equation; 

f = ^rT + 0,3 (1_5) 

Some of the most recent work in measurements of pressure drop across 





*Ref . 4, Fig. 217, p. 2-12. 



woven screens has been performed by Martin-Marietta of Denver (Ref. 6) 
using four dutch twill screens: 

• 375 x 2300 

• 325 x 2300 

• 250 x 1370 

• 200 x 1400 

with GN2 as a test fluid. Their resulting data, as shown in Figure (1-4)*, 
lie above the Armour-Cannon correlation for Reynolds number less than 10” 2 
and between the Armour-Cannon correlation and the MDAC test data for Reynolds 
number in the range 10” 2 < Re < 1.0. 

Information in the paragraphs above is summarized in Table 1-1 and 
Figure (1-5). 


*Ref. 6, Figure 11-38, p. 11-60 


Equations (1-2) and (1-4) were developed by GDCD (Ref. 7) after 


examination of their data arranged in the form f = a/Re + 3. They noticed 
that considerable error could arise in attempting to use a single 
correlation for all screens. GDCD, therefore, proposed that the most 
accurate way to arrange the data was as follows: 


A L y A L p 

Ap = v + V 2 


^a gc 


Da gc 


(1-6) 


where 


Ap is the pressure drop across the screen 
y is the fluid molecular viscosity 
p is the fluid density 

g c is the conversion constant needed if the "American engineering 

system" is used ( = 32.174 — ^~z ) 

\ lbf sec / 


V is the screen approach velocity 

D a , L, A]_, A2 are as defined in equation (1-4). 

Equation (1- 6) is merely a rearrangement of Equation (1—4) using the 
definition Eu = 2g c Ap/pV 2 and Re a = pVD a /y. Equation (1-6) may be further 
simplified in the following manner: 


Ap = A y V + B c p V 2 


AtL 

where A c = — - 


D a Be 


and E c = 


a 2 l * 
Da&c 


(1-7) 


GDCD presents values of A c , B c , A^, and A2 in Table 2-2 of Reference 
7. Values of the solidity S for square weave screen are presented in 
Table 2-3 of the same reference. It must be pointed out, however, that there 


*A C and B c are referred to as A and B in Reference -7. 


is a discrepancy in this data, as either or A 2 are not dimensionless or 
A c and B c are not in the units given in their paper. Tables 1-2 and 1-3* 
summarizes the values of A c , B c , A^, A 2 and S. 

North American Rockwell (Ref. 8 ) presented data for eleven screen 
materials, in the form given for equation (1-4) using liquid heptane as 
the test fluid. Over the range of Reynolds number tested, though, they 
found that their data could best, be arranged by deleting the second term 
containing A 2 inequation (1-4). Letting K= 2LA^/D a , the equation Eu = K/Re 
is obtained. Values of A and K are presented in Figure 1-6 

1-2 Experimental Equipment and Procedure 

Experiments were conducted on 200 x 200 and 325 x 325 square weave 
screens and 50 x 250, 200 x 600, and 325 x 2300 dutch twill screens using 
tap water as the test fluid. A schematic drawing of the first experimental 
apparatus is shown in Figure 1-7, and the screen/ channel assembly in 
Figures 1-8 , 1-9 and 1-10 . 

The experimental assembly consisted of the following pieces of 
equipment : 

(1) A centrifugal pump capable of producing outputs from 40 cc/sec 
to over 400 cc/sec. 

(2) A filter system capable of 5 micron filtration. This was 
used to remove impurities from the tap water and to act as a "buffer" to 

dampen the pulsating output characteristic of a centrifugal pump. 

*Ref . 7, Table 2-2 and 2-3, p. 2il4. 

+Ref. 8, Figure 4.5, 4.4-2, p. 
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(3) A rotameter, as shown in Figure 1-11, calibrated with tap 
water over the range of 30 cc/sec to 150 cc/sec. 

(4) A screen/channel assembly as shown in Figure 1-9, in which 
the screens listed above were set in place at location A. The channel 
itself was a standard 10.6 cm x 4.2 cm (4" x 1.647") channel size. A 
plexiglass tank, B, which was divided into two sections, sat on top of 
the channel, C , and inside screen mounting brackets, D, which were 
attached to the outside of the channel section (lower mounting bracket) 
and to the outside of the plexiglass tank (upper mounting bracket) . 

Water flow was introduced into one side of the plexiglass tank 

which contained an overflow tube, E, and then overflowed into the second 
section of the tank, which was directly over the screen. The second 
section of the tank contained a flow straightener , F, parallel to the 
screen. The tank thus served two purposes — the first section acted as 
a calming region for the inlet flow while maintaining a constant liquid 
head by means of the overflow tube, and the second section directed the 
flow perpendicular to the screen. Four sets of pressure taps were located 
in the screen mounting brackets directly above and below the screen. 

(5) A manifold system, shown in Figure 1-12, consisting of four 
three-way stopcocks connected to a single outlet tube. With this system, it 
was possible to close off three sets of pressure taps, and by rotating the 
remaining, open stopcock between its two openings, to measure the pressure 
above and below the screen. In a similar manner (after closing the first 
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stopcock), each of the other stopcocks could be opened, one at a time, 
and pressure differentials could be measured at each of the other 
locations quickly and with sustained accuracy. 

(6) An inclined tube manometer, shown in Figure 1-13, with a 0 
to 4 inch scale was used. The test fluid itself was used as the manometer 
fluid. 

(7) Other pieces of equipment (as shown in Fig. 1-7), including 

a water tank to furnish a liquid head for the pump, valves to regulate out- 
put and control liquid height in the tank above screen, and a 0° to 40°C 
thermometer graduated in tenths to measure the temperature of the water at 
steady state. 

The procedure for an experimental measurement was as follows: A 
volumetric flow rate was chosen by adjusting the pump outlet control valve 
(see Figure 1-7). The pump was allowed to run at this output until steady 
state was achieved . Steady state was indicated by the constant temperature 
of the water, constant output of the pump as indicated by the rotameter, 
constant liquid head above the screen as measured by a ruler fixed to the 
side of the tank (this height was the same for all screen tested) and no 
fluxion (i.e., the same reading) in pressure drop across the screen. 

Once steady state was achieved, the pressure differential at the four 
locations was recorded. A new, lower volumetric flow rate was then chosen, 
and the above procedure repeated. 

The upper volumetric flow rate was limited by the amount of flow 
a given screen would pass, while maintaining a constant fixed height above 
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the screen with the outlet control valve wide open. Increasing the flow 
rate beyond that limit merely increased the liquid head above the screen 
to such a point that it would eventually overflow the tank. 

The lower volumetric flow limit was fixed at the point at which 
the smallest scale pressure differential could be read (0.254 mm of water). 



The data collected during the test runs for the five screens are 

presented in Figures 1-14 through 1-17. It was found that the pressure drop 

Ap, was proportional to the average entering velocity, V , in the present 

e 

flow configuration (rectangular). Flow fields of previous investigators 

were in straight tubes or channels with screens pendicular to the flow 

*Note that the present work at UAH uses V e , which denotes a entering velocity, 
while previous investigations use V, which denotes a approach velocity. Thus V £ 
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>1 0 ) 


direction. 


Correlations of the present data against previous investigations 
have been made. 

Figure 1-18 is a comparison of the present data for square weave 
screens with the Hoerner equation, Eu =j^^J . Value? of the solidity, S, 
were taken as 0.70 and 0.66 for the 325 x 325 and 200 x 200 mesh screens, 
respectively. The equation Eu = 5.44 represents the 325 x 325 mesh screen 
and Eu = 3.77 represents the 200 x 200 mesh screen. The experimental data, 
for both screens, lie below the Hoerner equation- at high velocities and 
above the equation at the lower velocities. For the lower velocities, the 
slope of the Hoerner equation and that of the present data is nearly 
identical. 

A comparison of the present data, for the three dutch twill screens 
tested, against equation (1-4), was not possibly due to the discrepancy 
(as mentioned previously) of the units of A c and B c or the values of Ai and 
A'2 (given in Table 1-2 of this report) . A comparison of the data for the 
50 x 250 dutch twill screen with the North American Rockwell value of A^ 
was not possible because they tested a plain dutch screen. 

Present data (with the entering velocity changed to approach velocity 
for comparison with previous investigators) for the 50 x 250 and 325 x 2300 
dutch twill screens are shown compared to works by previous investigators 
over the same Reynolds number range in Figure 1-19. No value of "a" 
was available in the literature for the 200 x 600 screen. Screen parameters 
such as £, B, etc., available in the literature vary from source to source - 
those used in all calculations made for this report are given in Table 1-4. 
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The presently obtained plots are lower than any previously reported 
and are not linear. Two possible reasons are: 

(1) Geometric; i.e*, the present experimental apparatus is rectangular 
while previous investigations were "straight pipe flow" configurations. The 
Reynolds number as' defined by Armour-Cannon and GDCD does not take into 
account the channel diameter (or hydraulic diameter in the case of a non- 
circular configuration) but rather a characteristic screen pore diameter D a . 
Thus, Re, as defined, is independent of the type of geometric configuration, 
but it is doubtful that this is actually true. 

(2) Previous data were obtained from gas flow analysis while the present 
experiment used water as the test fluid. For a gas such as G^, the kinematic 
viscosity is an order of magnitude larger than water, thereby shifting the 
Reynolds number of water to the right (i.e. (Re) H 2 O > (Re) GN 2 ). Also for a 
gas, the friction factor would be less than that for water for the same Reynolds 
number. Thus, the present data viewed in light of previous gas data would 
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shift up and to the right, the exact magnitude being impossible to predict 
without simultaneous experimental comparison between gases and liquids. 

The present data for the 50 x 250 and 325 x 2300 dutch twill screens 
present a different tendency from former investigations - an increase in the 
friction factor with increasing Reynolds number up to a certain Reynolds 
number for both screens. At low velocities, the pressure drop for a fluid 
such as water is extremely low; thus, the relationship between Ap and V 
causes the friction factor to rise for low V (f* Ap/V 2 ) and fall for high 
V. It can be considered that previous data for gases would show this 
tendency if velocities were low enough to obtain the small pressure differences 
measured in the present case. If the data are analyzed only in the region 
where f decreases for increasing Re, equation (1-1) satisfactorily aligns the present 
data points for the two dutch twill screens. 

1-4 Recommendations for Future Work 

Examination of the experimental results suggests two additional areas 
of study: 

(1) Because of the geometric configuration of the present experimental 
apparatus, it is impossible for a given screen to obtain data over a wide 
Reynolds number range, as has been done in the previous "straight pipe-flow" 
experiments. With the exception of a very limited Reynolds number range, 
this prevents comparison of liquid data with a correlation of the form 
f = a/Re + 3 , which can be applied to a gas flow. It is required, 
therefore, that fluids such as water and heptane should be used in "straight 
pipe-flow” apparatus with the dutch twill screens tested in this experiment 
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in order to compare the previous gas data with measured values of liquid. 

The result would be the extension of the Reynolds number range of liquids 

over that of formerly tested gases. It would allow comparison of the present data with 

MDAC (Ref. 3) and GDCD (Ref. 4) test data, and, finally, a comparison could 

\ 

be obtained for the same screens in different^geometric configurations 
(i.e., rectangular and circular). 

(2) For velocities, V e , used in the present experiment (approximately 
1 -* 5 cm/sec), each dutch twill screen has a Reynolds number range of less 
than 10. New dutch twill screens could be chosen to cover Reynolds number 
range different from those tested in the present case [50 x 25Q Re 0.6 -* 3, 

325 x 2300, Re 0.03 •* 0.06]. This would accomplish three objectives: 

• It would allow the "holes" in Reynolds number to be covered, 
thus providing a correlation of liquids which would be 
indicative of the whole Reynolds number range anticipated. 

• It would furnish additional or new data on many dutch twill screens. 

• It would allow a more complete correlation with previously collected 
gas data. 

In addition, more experience would be gained in the improved design and use of 
rectangular-channel screen assemblies. 
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SECTION 2 

VARIATION OF PRESSURE DROP ACROSS WOVEN SCREEN 
2-1 Theoretical Analysis of Effect of Nonunlform Flow, McDonnell Douglas 

Previous investigators have developed correlations to predict pressure 
drop across woven screens as a function of a uniform and average approach 
velocity (volumetric flow rate divided by screen area). However, the division 
of a fluid stream into parts by means of a screen is accompanied by pressure 
changes owing to friction and the change of fluid momentum. As a result, the 
pressure drop is not constant along" a screen and may be considerably higher 
than that predicted by an equation based on an average approach velocity. 

McDonnell Douglas (Ref. 9, Appendix A) considered the channel 
configuration shown in Figure 2-1.* The following equations were presented; 

Continuity 
where 

Momentum 

where 


*Ref. 9, Fig. A-i, p. 92. 


d u 

V = D (2-1) 

dx 

V is the velocity normal to the screen 
D is the channel depth 
u is the velocity in the x-direction 

g +f _£_ ^ !p ,| . ° (2-2) 

p is the static pressure 

f is the friction factor defined by Darcy's equation 
Djj is the hydraulic radius 

g is the acceleration due to gravity 
p is the fluid density. 
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The flow loss through the screen was assumed to be 

KqV - P 0 - Pgx - P (2-3) 

where K 0 is determined experimentally. Equation (2-3) is applicable for low 
velocity and in terms of the Armour-Cannon correlation is valid where 8. 6/Re » 
0.52 (g, an inertial resistance coefficient, defined in equation (L-l), is 
negligible) . 


With the following boundary conditions and definitions, the three 
equations above are non-dimensionalized 


B.C. 

(1) 

at x = -H 

u = 

0 

(2-4) 

B.C. 

(2) 

at x = 0 

Q = 

pLDu 

(2-5) 

H is 

the 

channel length 




L is 

the 

channel width 




Q is 

the 

mass flow rate 







u 

PuDL 

(2-6) 



u* — 

Umax 

Q 



V* a 

V 

pVLH 

(2-7) 



V 

avg 

Q 





Z = 

x + H 
H 


(2-8) 


Ap* = 


Po-P 
Ap 


avg 


(Po-P) PLH 

K 0 Q 


(2-9) 


where ■ the terms which contain asterisks are dimensionless and Z is the 
dimensionless length. 

The three basic equations (2-1, 2-2, 2-3) become 
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v* 


du* 
“d T 


(2-10) 


Ap* - 


p 2 LH 2 

"k^q" 


£ (z-1) 


d(Ap*) fQH 2 , p 2 LH 2 e 2QH du* 

u * 2 - “ — = , U* — ; — 

dZ 2DhD 2 LKo KqQ D 2 LK q dz 


( 2 - 11 ) 


which are combined Into a single nonlinear equation 

^ - F 0 u* 2 - 0 u* ^ = 0 (2-12) 

where F = fH/4D h and 0 - 2QH/D 2 LKo- 

The boundary conditions are written as follows: 

B.C.(l) at Z = 0 u* = 0 (2-13) 

B.C. (2) at Z = 1 u* = 1 (2-14) 

Equation (2-12), with the assumption that F = 0, is solved by McDonnell 

Douglas. The results are presented in Figure 2-2* as a plot of V* versus 

0. The figure shows that V* at the ends of the channel begins to differ 

significantly from 1.0 at values of 0 greater than 1.0. In other words, V 

is not the average approach velocity for 0 > 1 (for <j><l, V* = r~ — ~ 1.0). 

v avg 

2-2 Extension of McDonnell Douglas Analysis of Effect of Nonuniform Flo w 

McDonnell Douglas solved equation (2-12) with the assumption that F « 0. 
It was felt that the solution for equation (2-12) would be different by 
including a treatment of F, especially In the case of small channel depth, 

D, or high volumetric flow rate Q' (corresponding to high velocities V). 

*Ref . 9, Fig. A-3, p. 98. 
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A computer program employing an implicit finite difference technique 
of Crank-Nicholson type was developed to solve equation (2-12) in this 
study. The computer program listing is included in Appendix B. The results 
of the calculations are shown in Figures 2-3 and 2-4 as plots of V* versus 
0 at Z - 0 and Z = 1 with the friction factor in the x-directisn, F , as a 
parameter. Other plots of V* versus 0 and V* versus Z with F as an independent 
parameter are presented in Appendix B. 

A comparison of Figure 2-2 with Figures 2-3 and 2-4 reveals the fact 
that the effect of F is negligible in the range F _< 1.0 and may be neglected 
as originally proposed by McDonnell Douglas. 

2-3 Experimental Apparatus and Procedure 

The experimental apparatus previously discussed in Section 1-2 and 
illustrated in Figures 1-7 through 1-13 was modified to investigate the 
effect of nonuniform flow (nonuniform flow means that the local entering 
velocity, V e .varies along the screen because of the accelerating effect of 
the velocity u which is parallel to the screen) . The following changes 
were made: 

(1) The channel length H.was increased from the initial length of 
8.89 cm to a new length of 29.16 cm. 

(2) The channel depth was variable by the insertion of plexiglass 
spacers which reduced the depths of the channel from an original depth of 
3.5 cm to new depths of 2.85 cm and 2.2 cm. 

(3) A total of fourteen sets of pressure taps were located at the 
following positions where the origin x = 0 corresponds to the end of the 
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channel as shown in Figure 2-5. 


Location Distance x, cm Dimensionless Distance Z *» 


1 

0.28 

0.0096 

2 

1.79 

0.06 

3 

5.95 

0.20 

4 

9.08 

0.31 

5 

10.69 

0.37 

6 

12.32 

0.42 

7 

13.97 

0.48 

8 

15.53 

0.53 

9 

19.56 

0.67 

10 

21.17 

0.73 

11 

22.76 

0.78 

12 

24.37 

0.84 

13 

27.18 

0.93 

14 

28.67 

0.98 


(4) A new rotameter calibrated with water over a range of 100 cc/sec 
to 400 cc/sec was installed to allow higher volumetric flow rates . This 
resulted in obtaining the same average entering velocity range, V e , in the 
modified set up as in the original experimental apparatus. 

(5) A larger outlet was provided at the end of the channel to 
accommodate the higher volumetric flow rates. 

All other pieces of equipment shown schematically in Figure 1-7 were 
retained. A photograph of the modified screen/channel assembly appears in 
Figure 2-6. 
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The experimental procedure used in the modified set-up was the same 
as described in section 1-2 with the following exceptions: 

(1) Pressure drop measurements across the screen were taken over a 
volumetric flow rate range of 100 cc/sec to 400 cc/sec instead of 30 cc/sec to 
150 cc/sec. 

(2) Each volumetric flow rate chosen was repeated four times in 
order to measure the pressure drop at the fourteen locations (only four 
locations at a time could be measured with the manifold system shown in 
Figure 1-12). 

2-4 Experimental Results 

Data were taken for a 50 x 250 dutch twill screen using tap water in 
the modified experimental apparatus discussed in Section II-3. The volumetric 
flow rate, Q^, was varied from 400 cc/sec to 100 cc/sec. The depth of the 
channel, D, was changed from 3.5 cm to new depths of 2.85 cm and 2.2 cm. The 
pressure drop across the screen was measured at the fourteen locations. 

Figure 2-7 presents the pressure drop as a function of volumetric 
flow rate at the two ends of the channel (Z = 0.01 and Z - 0.98) for depths 
of 3,5 cm, 2.85 cm, and 2.2 cm. The effect of the channel depth, D, on the 
pressure drop, Ap, is considerable. At Z = 0.98 the channel depth plays a 
significant role in increasing the pressure drop for high volumetric flow 
rates (which induces the higher velocity u) . At Z = 0.01 where the velocity 
u is almost negligible, the smaller channel depth reduces the pressure drop 
in contrast to the result at Z “ 0.98. 

Figure 2-8 presents the same pressure drop data, as in Figure 2-7, 
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plotted against the average entering velocity V . The average entering 
velocity is calculated from V g = Q'/A& where e is the screen void fraction. 
Since the average entering velocity and the local entering velocity are not 
identical. Figure 2-7 is a more accurate representation of the data than 
Figure 2-8. 

Figures 2-9 through 2-21 present the pressure drop as a function of 
the dimensionless length, Z, for different volumetric flow rates at depths of 
3.5 cm and 2.2 cm. 

Analysis of the data of Figures 2-7 through 2-21 yield the following: 

(1) The average entering velocity, and the local entering 
velocity, V e , are considerably different at the higher volumetric flow 
rates (Q '> 200cm 3 /sec) and extreme end (Z = 0.98) of the channel. In terms 

of the McDonnell Douglas correlation, this means that both V* and 0 are greater 
than one . * 

(2) The relationship between Ap and V e is no longer linear over the 
same average entering velocity range as studied in phase one of the experi- 
mental work (Section 1-3). 

(3) Pressure drop across the screen is a strong function of channel 
depth D. At the smallest depth of 2.2 cm, the pressure drop difference be- 
tween the two extreme ends of the channel (Z = 0.01 and 0.98) is the greatest 
while at the largest depth of 3.5 cm, the difference between the two ends is 
the smallest. If the channel depth was increased to some finite value, 

the pressure drop difference between the ends of the channel would go to zero 
and then V e and V e would become identical. 
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2-5 Analysis of Experimental Data 


2-5 . 1 McDonnell Douglas Correlation of V* versus 0 

The experimental data of Figure 2-8 were first plotted in the form 
of V* versus 0 in order to analyze the McDonnell Douglas theoretical model 
discussed in Section II-l. Results are presented in Figure 2-22; however, 
the parameters V* and 0 are modified in the following manner: 

(1) V is based upon a fictitious velocity, not a local entering 
velocity; that is 


V = p(V e ) at Z LH/Q 
atZ 


where 


(V e )at Z = 


AP at Z 
(Ko) av g 


(2-15) 

(2-16) 


(2) 0 is based upon an average value of Kq; that is 


0 - 


2QH 

D 2 L(K 0 ). 


where (K 0 ) avg 


o ' avg 

( K o)at Z = 0 + < K o)at Z 


= 1 


(2-17) 

(2-18) 


Equation £-18) is used in place of the original McDonnell Douglas definition. 



(2-19) 


Equation (2-19) implies a linear relationship between pressure and velocity; 
an inspection of Figure 2-8 shows that this is not the case. Therefore, 
equation (2-16) is used as the definition of K 0 ; that is, K 0 is Ap at Z divided 
by V e at Z and since Ap at Z varies as a function of Z, then Ko is a variable 
and a function of Ap at Z , V g at Z , and Z . 

V* and 0 in Figure 2-12 are calculated in the following manner. 
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(1) For a given volumetric flow rate, an average entering velocity 


is calculated. 

(2) For the average entering velocity calculated in (1) , Ap is read 
from Figure 2-8 at Z = 0.01 and Z = 0.98 and the value of K 0 are calculated 
at Z = 0.01 and Z = 0.98. Than (K 0 ) aV g is calculated by Equation (2-18). 

(3) Once (Ko) av g is determined, V e is calculated at Z = 0 and 1 
(actually Z = 0.01 and 0.98) by equation (2-16). Then V* is calculated at 
Z = 0 and 1 by equation (2-15) . 

(4) 0 is calculated by Equation (2-17) using (Ko) aV g determined in 
step (2) above. 

Thus, a given value of 0 yield two values of V*, one at Z = 0 and Z = 1 

The results shown in Figure 2-22 follow- the same trend as the McDonnell 
Douglas theoretical curve. However, the divergence of the curves for Z = 0 
and Z = 1 at a given 0 is not as great as predicted by theory. Finally, the 
plot in Figure 2-22 is based upon entering velocities not approach 
velocities; thus the curves converge to a value of approximately three 
instead of one. Since V e = V/e and the void fraction for a 50 x 250 dutch 
twill screen is 0.325, Figure 2-22, if based upon approach velocity, would 
shift down to a value of approximately V* = 1. 

Several points should be made about Figure 2-22, the calculation 
method used to determine V* and 0, and the McDonnell Douglas theoretical 
model in general. 

(1) The results of Figure 2-22 are* encouraging in that the simple 
equation Ap = K 0 V may be used to predict the variation in pressure drop 
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as a function of Z. This is true because in the present measurement the 
local static pressure above the screen was almost uniform for a given 
volumetric flow rate and channel depth while the local static pressure 
below the screen varied in the Z-direction. This variation caused the 
pressure drop difference. 

(2) The validity of the equations used is questioned, A linear 
relationship is used for Ap even when the actual relationship is quadratic. 
Also, a great deal of ambiguity exists in both the definition of Ko and the 
determination of its value. Finally, without recourse to an actual measured 
or calculated velocity profile, a fictitious velocity profile must be 
calculated. 

(3) Realizing that the definitions employed for V* and 0 are not 
those originally proposed by McDonnell Douglas we still feel that the data 
as presented in Figure 2-2 or Figure 2-22 are not a good way to predict the 
variations in pressure drop. Two reasons are proposed. First, the correlation 
is based upon the fact that the inertial contribution to the pressure drop is 
negligible. With this assumption, a linear relationship Ap = KoV + pgx 

(or Ap = K 0 ) is defined - this is perfectly valid since the viscous 
contribution is assumed to be the only contribution to the pressure drop. 
However, the present experimental data for the screens tested has shown that 
when a linear relationship exists between Ap and V, no discernible pressure 
difference can be measured. Thus the local velocity and average velocity 
are identical. As the pressure drop difference at Z = 0 and Z = 1 begins 
to differ significantly, the assumption that the inertial term is negligible 
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is no longer valid and the defining relationship for K 0 should include a V 2 
term. Second, even in the viscous range, when the defining relationship 
for K 0 is valid, the question shall arise as to what value of Kq to use - 
a great deal of ambiguity exists for this term. 


2-5.2 Determination of Velocity Profile 

The equations of continuity and momentum are derived, in detail, for 
the present experimental apparatus in Appendix A. The equations are: 


Continuity 


Momentum 


Ve = 


D du 


e dx 


lE + f£Hl + 3pu^i = 0 

dx 2D dx 


( 2 - 21 ) 


( 2 - 22 ) 


Assuming that the second term of equation (2-22) can be neglected as noted in 
section 2.2, then equation 6-22) is written as 


dp du n 

+ 3 pu — = 0 

dx dx 


(2-23) 


Integrating equation(2-23) from zero to x 


u(X) = ^ 


2[Ap(x) - Ap(0)3 
3 P 


(2-24) 


where at x = H, u(H) = — 


pLD (2-25) 

Now, using equation (2-24) and experiment data at X = H, u(H) is calculated. 

The values of u(H) calculated by equation (2-24) are compared with the values 

calculated by equation (2-25) which is exact. A correction is made to 

equation (2-24) so that the values calculated at x ■ H agree closely with 

those calculated by equation (2-25) . Thus equation (2-24) becomes 
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(2-26) 


ufc) = 1.18 

and is used for all subsequent calculations of the velocity u at location x. 

The velocity* V e , is now calculated from the continuity equation once 
du/dx has been determined from Figure 2-23. 

Figures 2-23 and 2-24 show the velocity profiles u(x) and V £ versus 
Z for a 50 x 250 dutch twill screen in a channel of depth 2.2 cm. 

By comparing V e versus Z (Figure 2-24) and Ap versus Z (Figures 2-16 
through 2-21) for D = 2.2 cm, Ap as a function of V e is obtained. This is 
shown in Figure 2-25. The data, for low velocity, is nearly independent of 
the mass flow rate Q. 

For practical calculations, a single line can be drawn through the 
data points as shown. 

Based upon the single line of Figure 2-25, Euler number is plotted as 
a function of Reynolds Number in Figure 2-26. An empirical relationship between 
Eu and Re may be obtained from this figure. 

2-5.3 . Determination of Average Properties from Local Data for 50 x 250 
Dutch Twill Screen for Depths of 3.5 and 2,2 cm. 

As most existing correlations are based upon the assumption of an 
average and uniform velocity, the present data were analyzed to determine the 
average pressure drop, Ap, as a function of the average entering velocity, V e . 

Figures 2-9 through 2-15 were integrated at each volumetric flow 
rate to determine Ap. Each volumetric flow rate corresponded to an average 
entering velocity V e (V e = Q /Ac). Thus, Figures 2-27 and 2-28 present Ap 
versus V e for depths of 3.5 cm and 2,2 cm, respectively. 
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Figure 2-29 presents the Euler number as a function of the Reynolds 
number based upon the average properties of Figures 2-27 and 2-28. At low 
Reynolds number, the Euler number is a function of the depth of the channel. 

As the Reynolds number increases, Eu becomes independent of channel depth. 

Other dimensionless parameters, such as those shown in Figure 2-30, 
were investigated to see if the data for different channel depths, D, would 
conveniently clasp into one single curve. Figure 2-30 was found to be 
the best possible presentation which includes the depth, D, and channel length, 
H, as parameters. 

2-6 Recommendations for Future Work 

Based upon an analysis of the data of Section II , the following 
points are made and suggested for future study. 

(1) Since only one available screen (50 x 250 dutch twill) and 
one liquid (water) were tested, it is obvious that other screens and 
liquids should be tested in the experimental apparatus discussed in 
Section II-3. Special attention should be paid to the viscous region 
where it appears that such parameters as Euler number are dependent upon 
channel depth. 

(2) It is recommended that an approach to pressure drop prediction 
as discussed in Section II-l be taken only as an approximation as any 
single correlation for all screens can not be accurate enough for detailed 
calculations. It is believed that the following equation is the most 
accurate way to represent the data for each single screen. 
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2L 

Al 


Eu 

ss ' 

— + A* 

(2-27) 


Da 

Re 

or in alternate form Ap 

= a c p v + 8c y v z 

(2-28) 


If further studies verify the fact that Eu is a function of channel depth in the 
viscous region, then equations (2-27) or (2-28) must either present value of A c 
and B c (or and A 2 ) as a function of depth or modify the definition of 
Reynolds number to include the effect. Based on the present work, it 
appears that equations (2-2 7) and (2-28)may be used in the inertial region 
without modification. 

(3) If Equations (2-27) and (2-28) are to be used, an accurate 
determination of a velocity profile is needed. Accurate velocity profile 
should be measured for the screen tested (50 x 250) and others. The 
measured velocity profile for a 50 x 250 dutch twill screen should be 
compared with the calculated velocity profile (as described in II-5.2); 
if the two are in close agreement, then velocity profiles may be calculated 
from the following; 


u (x ) = n 


Apfr) - Ap(0) 

P 


where n is a constant 


V e 


I) du 
e dx 


(2-29) 


(2-30) 
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TABLE 1-1: SUMMARY OF RESULTS OF INVESTIGATORS 

USING ARMOUR-CANNON TYPE CORRELATION 



Investigators 

or 

Company 

Equation Developed 
or Results 

Screens Tested 

Fluid Tested 

*3 if ■ 

Armour-Cannon 

f - 8,61 + 0.52 

Too Numerous to 

GN 0 , 

— 

(Ref. 2) 

RS 

List 

GHe 


MDAC 

No equation given; 

■ 



(Ref. 3) 

results below 

250 x 1370 


u a 


Armour-Cannon 


gn 2 . 



generally by a 

325 x 2300 

GHe 



factor of two. 



j * 


Armour-Cannon 

200 x 1400 




correlation sue- 





cessful in aligning 




..... J 

data points 




GDCD 

• 

f - 2>/+9 + 0.3 

Six screens 

gh 2 


(Ref. 4) 

Ke 

tested in 1969, 

gn 2 

w * 



Three screens 

Too numerous 

■’ 



tested in 1968* 

to list 

* M 

Martin- 

No equation given; 

375 x 2300 



Marietta 

results above 

325 x 2300 



(Ref. 6) 

Armour-Cannon 

250 x 3700 

gn 2 



correlation for 

200 x 1400 




Re < 10 -2 results 



IU 


between Armour- 





Cannon correlation 





and MDAC data for 





10- 2 < Re < 1. 
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Table 1-2*. Screen Geometry, Dutch Twill Screen 


Screen Mesh 

°B 

Bubble Pt. 

Diameter 

(Microns) 

»a d * 

Avg. Capil- 
lary Dia. 
(Microns) 

Thick- 

ness 

(in.) 

Por” 
osity 
(Meas . ) 

Dimensionless 

A.L 

A — * 

A 2 L 

H 

1 

Wicking 

A w 

A o 

C z 

Dg 

a e c 

B c D g 
a c 

200 x 1400 

13.4 

I 

0.0058 

0.256 

llll 

18 


509687 

3.61 

165 x 1400 

18.6 


0.0060 

0.301 


16 

580 

270185 

2.68 

200 x 600 

19.05 

36.7* 

0.0055 

0.347 

52 

3 

368 

51053 

0.355 

165 x 800 

22.7 

48.5** 

0.0065 

0.310 

43 

135 


28568 

14.3 

150 x 700 

22.7 

60.8** 

0.0070 

0.171 

S|m 

133 


227642 

12.08 

80 x 700 

29.7 

139.3** 

0.0098 

0.416 


34 

6230 

121427 

1.89 

50 x 250 

33.9 

129.5** 

0.0127 

0.325 

115 

191 


20938 

14.78 

30 x 250 

48.5 

112.2** 

0.0265 

0.276 

130 

12 

1120 

65795 

2.23 

20 x 250 

52.8 

155.3** 

0.0280 

0.325 

150 

20 


41869 

2.84 

♦Microporosimeter 









♦♦Macroporosimeter 

i 









Table 1-3: Screen Geometry, Square Weave Screens 

•v 


Mesh 

(Microns) 

d b 

Solidity 

S 

400 x 400 

• 38 


325 x 325 

44 

1 

200 x 200 

74 

0.66 

150 x 150 

. 104 

0.63 

100 x 100 

140. 


80 x 80 

180 

1 

50 x 50 

280 

0.700 

40 x 40 

440 

0.640 

20 x 20 

860 

0.538 
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TABLE 1-4: GEOMETRIC PARAMETERS FOR WOVEN SCREENS 


Screen 

D a , cm 

B, cm 

a, cm -1 

e 

B 

S 

200 x 200 

- 

- 

- 

- 

- 

0.66 

325 x 325 

- 

- 

- 

- 

- 

0.70 

50 x 250 

12.95 x 10- 3 

32.3 x 10 -2 

151 cm" 1 

0.325 

1.3 

- 

200 x 600 

3.67 x 10 -3 

14 x 10- 3 

- 

0.347 

1.3 

- 

325 x 2300 

1.472 x 10- 3 

8.9 x 10- 3 

1102.3 

0.245 

1.3 

- 
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The table below summarizes the changes in the experimental apparatus 


TABLE II- l: MODIFIED EXPERIMENTAL APPARATUS 


Piece of Equipment 

Modifications 

Screen/Channel Assembly 

(1) length H increased from 8.89 cm 
to 29.16 cm. 

(2) width L unchanged (8.89 cm) 

(3) depth D variable (3.5 cm, 

2.85 cm, 2.2 cm) 

(4) number of pressure tap 
locations increased from four 
to fourteen 

(5) outlet increased in size 

centrifugal pump 
filter 

inclined tube manometer 
control valves 
feed tank 
manifold system 

none 

screens 

only 50 x 250 dutch twill screen 
used 
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Figure 1-2 Comparison of MDAC Data with 
Armour-Cannon Correlation 
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Figure 1—7. Schematic of Experimental Apparatus 




Figure 1-8. 
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Figure 1-12. 
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Figure 1-13. 
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Figure 1-14 Pressure Drop as a Function of Average 
Entering Velocity, Square Weave Screens 
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Figure 2-6. Modified Screen/Channel Assembly. 
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Figure 2-8 
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Figure 2-10 

AP versus Z, Q = 350 cm 3 /sec, D = 3.5 cm 
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Figure 2-13 

AP versus Z, Q - 200 cm 3 /sec, D - 3.5 cm 
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Figure 2-14 AP versus Z, Q - 150 cm 3 /sec, D - 3.5 cm. 
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Figure 2-15 AP versus Z, q = 100 cm 3 /sec, D - 3.5 cm 
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Figure 2-18 AP versus Z, Q - 300 cm 3 /sec, D - 2.2 cm 
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Figure 2-21 AP versus Z, Q - 100 cm 3 /sec, D - 2.2 cm 
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Figure 2-25 Pressure Drop as a Function of Entering Screen Velocity 
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Reynolds Number, Re = — *— 

a v 


Figure 2-26 Euler Number versus Reynolds 
Number (Based on Figure 2-25) 
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Figure 2-29 

Euler Number Versus Reynolds Number (based on average pressure 
drop and average entering velocity). 
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APPENDIX A 


This appendix contains the complete derivation of the continuity 
and momentum equations for the experimental model used at UAH. 

The notation used in this appendix is the same as that used in 
Sections 1 and 2 of the report and listed in the definition of symbols. 
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m - (m + — Ax) LD - puVLAxe - p + ^ Ax-p LD+f'LAx (4) 
dx dx 


simplifying (4) 


" dx (Pu 2 )AxLD- puVLAxe = ^ Ax LD + f'L Ax 


(5) 


dividing both sides by AxLD 


d , 
— (Pu 2 ) 
Qa 


puVe 

D 


<&, L 

dx ^ D 


( 6 .) 


rearranging 


^ + 2pu— + 


uV 


dx 


dx 


D 


= 0 


(7) 


du Ve 

from equation (2) <jx ~ D 


d£ + 3pu du + r 
dx dx D 


= 0 


( 8 ) 


now, let f" = f 


pu_ 

2 


then 


d£ , fP u . 0 „ du 
+ — rr — + 3 pu , =0 

dx 2D dx 


(9) 




:2 



I ■ ■ i 1 


Let us define 
t' = f 


Pu 


then Equation (6) is rewritten as 
dP , du , pu 2 _ _ 

s + 3pu di * £ V - 0 


Nondimens ionalizing Equation (8) 


(7) 


( 8 ) 


Equation (8) is written in a nondimensionalized form as 


d 2 u* 


- F 0 u* 2 - 0 u * 


du* 


dz" " " " “ dz 

with the boundary condition 

u* = 0 at Z = 0 

u* = 1 at Z = 1 


= 0 


(13) 


(14) 

(15) 


Let us omit the asterisk, *, from now on for convenience. Equation (13) 
will be numerically solved by utilizing the implicit finite difference 
technique of Crank-Nicolson as briefly introduced in the next section. 


Numerical Solution 

Defining the stations Z^, and Z it ^ as n-1, n, and n+1, 

respectively, as shown in the figure below: 
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we can write the first and second derivatives of velocity as 




du ^n+1 - ^n-1 


dZ 

d U ~n+l 


2AZ 

U - 2U + U 


n n-1 
1 


AZ 

Thus Equation (13) is expressed in a finite difference form as 

U . ' - 2U + U n U n+1 - U n „-, 

- F0 U 0' - 0D' -^rz = 0 

AZ^ n n n 2AZ 

where is an old value of U n obtained by a previous iteration at 

station n. 


( 16 ) 


(17) 


(18) 


)I|, 


We solve the equation below: 

Al U n -1 + A 2 U n + a 3 U n+ i = D n 

where 

Coefficients of U n +i: 
i 

A 3 ~Kz T ~ 2AZ 


(19) 


( 20 ) 


i 


1 

:* ; 


III 


Coefficients of U n : 


A 2 ‘W ' "n 


Coefficients of U r _^ ; 

0U 

r 

AZ* ' 2AZ 


a 1 , " n 

= -rzrr + 


D n = 0 


( 21 ) 


( 22 ) 

(23) 


Ml" 
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The boundary condition at Z = 0 is = 0, thus Equation (19) gives 

A(1,1)*U(1, JN)+A(1, 2) *U(2, JN) 

4A(1,3)*U(3, JN) = 0 
U (1, JN) = 0 

Therefore, we can set at Z = 0 
A(1 ,2) = A(l,2) 

A(l, 3) = A(l,3) 

A(l,l) = 0.0 


At A = 1, U n 


D(l) = 0.0 
= 1 , then 


max 


(24) 

(25) 


(26) 


A(NY2 , 1) *U (NY2 , JN)+A(NY2 , 2) *U (NT1JN) 

+A(NY2,3)*1.0 = 0 

Therefore 

A(NY2,1) = A(NY2,1) 

A(NY2 ,2) = A(NY2,2) 

D(NY2) = -A(NY2, 3) 

A(NY2, 3) = 0 

In the computer program the following symbols are defined; 
U(i,JN) = U n 
u(i,j0) = 

PHI -0 

i)z = Az 

VSTAR = dUn/dZ (=dU*/dZ) 


(27) 


(29) 
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For a given combination of constants, F and 0 the distribution of U, 
along Z is calculated. Iterations are continued until a desired 

convergence 

is obtained. The definition of convergence is 


u u' 

u n max u n max 



2 2 
U n max 

£ CONV. 



APPENDIX B 


CAPEFT Program Listing 

OFOR.lS MAIN. MAIN 

C0HMPN/BLUE/U<25-C,2» .OZ.NY.MVI ,Nv2 ,2 « ?5C I , D f 2B0 » . A I , 3 * 
COMMON /PTNK/PHl ,F 
COMNON/WH l YE/ JO.jN 

c 

N Y * 5 C 

NY I * N Y “ 1 

NY2*NY-2 

D P H I = £ • I C 
JN* J 
J 0 = 2 

D 7 * I ,C /FLOAT ( N Y 1 ) 

C 

1 F s C, . |(; 

2 PH I *0.10 
S CONTINUE 

DO If' ! * 2 • N Y 
I c Z(I)*07*FL0AT(I-I) 

2 ( I ) > l • 0 
C 

c U -PROFILE TO INITIATE CALCUtATlOM 

C 

DO 2 r, I = » * N Y 
?L U< I , JC>»Z< I ) 

C 

CALL &FEEL 
call print 

I F ( P H I . L T • » . 0 ) PH J =PH ! ♦DPP I 
I F ( P HI • G E • I * W ) PHi»PHI*l»C 
I F ( PH I • L E • i 5 #0 1 Go TO 5 
FsF + O. I 0 

I F ( F • GT • I • 0 ) STOP 

GO TO ? 

E M D 

ISFO R , I S GPEFN .GREEN 

SUBROUTINE GPEEN 

C COMMON/Rl UE/U ( 2SL , 2 > » OZ.NY ,NY I ,Nv?.Z<?SD) , D ( 2S0 > . A ( 2 5 0 t 3 i 

COMMON /PINK/ P H I , F 
COMMOM/WH I YE/ JO.jN 
M.NY/2 

NMAX = 2 SC 
I TMAX *■ I 00 
I T * t 

Z 1 S l • > / (OZ*OZ) 

: Z2 = 0.S0/DZ 

I pD CONTINUE 
la 2 ■ : ■ 

JO A ( I - 1 . 1) E ZI + Z?*PHI*UM,vJO) 

A * I - I , 2 V«-2.Q*Z I - pH I « U< ! . JO) *F*2»0/3»C t 
A C I - 1 ,3 ) «=* A ( I - I , I ',+7.0*7 I J 

0 (I - 1.1 PC -C 

1 = 1 + 1 



lr ( ! * L E . N Y l > GO TO 10 

C MODIFY FIRST AND (AST EQUATIONS BY B.C. 

A ( i , | ) * 0 » o 
D ( 1 ) ® • 0 

0 < N Y 2 ) = - A ( N Y 2 , 3 > 

A(NY? ,3) rC'.ii 

C 

C SOLVE EQUATIONS FnN U ( I , J '" ) , l * 2 • 3 • • • *' V I 

C 

Call trim ( A»U(2. JN) ,D»N v 2,NM<yX> 

c 

c applyR»c« for uri.JN) and u(ny,J m > 
c 

U( i , JN>*0*0 
U ( NY , JN ) *1 • 0 
C 

TEST* ( U ( f . JN ) -U ( M , JO > ) /U ( N i JN ) 

1 F ( ApS < TFST ) . LE .C. 0 I C.OR . f T . GT . I TMAX ) GO TO **C 

IT * IT 4 I 

DO 2 J K * I » N Y 
2 (i U ( K i JC ) *U ( K , JN) 

GO TO 100 
PD WRITE (&, 30) IT 

30 FORMAT ( 1 OX , 2>iHNlVBER OF I T ER A T I 0 N « r I B / ) 

RETURN 

END 

©FOR. IS TRIM, TRIM 

SUBROUTINE TRIM {-At X *0 • N , MN ) 

C 

DIMENSION A ( N N , 3 ) i A A ( 2 S r ) » 0 ( N N ) , D D < ? S C ) , X ( N N ) 

C 

C *• Vfppward TO ELIMINATION 

c 

A A < l ) * A ( I * 3 ) / A ( I . 2 > 

D D U ) * D ( I ) / A ( "t , 2 ) 

DO | I * 7 • M 
A A A « A t I , 2 > - A A ( I *- I 
AAU)=ACJ,3)/AAA 
C 

1 DD(M = (n{IJ-DO(I-t)*A(I,I))/AAA 
C 

c •*« pack Substitution • • * 

c 

x < n > = r> d < 

DO 2 I*? * N 
J»N- 1 * I 

2 X ( J ) * D D ( J > -X ( J + I > * A A I J > 
return 

END 

©FOR, IS PRINT, PRINT 

subroutine print 

COMMON/ BLUE /u ( 250, 2 ) , DZ , M Y , N Y t , N Y2 , Z < 250 ) , 0 ( 250 > , A (2 $ 
COMMON /P I Mk /P PI, F 
COMMON/WHl YE/ J0,JN 


WR I TE < A , t 0 » F , PH I 
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It FORMAT ( I H 1 , I L X , <Uip = , OFF I I . S , l 0V6HPH ! = ,OPEII.£>/) 

WP I TF < A , ?C ) 

pC FORMAKhH N 0 « , I ti X t| H 7 ,»f,X9M V E I OC 1 T Y , I £ X 8 H 0 L 0 V E L • , ! 0 X 6 H V S T A K ) 

1 = 1 

3C CONTINUE 

IF ( ! iGT, I . AMT . I .LT.NY ) VSTAR = U.«;0* I u U ♦ I , JN ) -U < I -1 , JN > ) /DZ 

I F ( t ,FG . 1 ) VST AP = ‘ . 50* ( - U ( 3 , JN' ) + « . *IJ < ? , JfM - 3 . *U I 1 , JKI ) ) / nZ 

l F ( I . FQ . MV ) VSTA^s < U ( MY , Jm ) - |i ( hjY | , JN M / ny. 

WR I TF. ( 6 , I Ci ) I , 7 ( I ) ,U ( I ,.JM ) ,H ( I ,,I0 ) , V-5TAR 
iro F 0 R F A T < I S * C; F R F I A • 7 ) 

I = I ♦ 1 

I F ( T *■ 6 T i H Y ) P F TOR N 

GO TO 3'C 

END 
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Sample Outputs 


f ■ .loooo phi ■ .loooctno 



Z 

VELOC ! TY 

OLD VEL • 

VStaR 


• OOOOqOn 

. ooooooo 

♦ ooooooo 

.90270 1 0 


•020400? 

.0200560 

.0200520 

.9827813 


» 0 <i 0 8 163 

.0001135 

#0001008 

.982341 9 


.06 1 2245 

.0601723 

• 0601590 

.9029030 


.0816327 

.0802336 

•0002160 

.9830809 


• 1 020HOP 

. 1002983 

• 100276R 

.9832676 


. 5 2 2 4 0 9 p 

. 1 203670 

. 1 2038 1 0 

•98309J5 


. 1 8 2 8 5 7 1 

. 1 0 0 0 0 C 8 

• 1 800 1 1 1 

.9837567 


.1632653 

.1605208 

• 1 600068 

.9800630 


# 18 3 6 7 35 

. 1 806066 

. 1005692 

#90001 1 B 


. 2n808 1 6 

. 2(307000 

• 2006593 

.9800022 


• 2200090 

.2208026 

•2207500 

.9052300 


. 2 0 0 8 9 0 D 

.28091 0 1 

#2008660 

.9857 100 


.265306 1 

.2610357 

#2609800 

.9862279 


* 2857 1 83 

.2811693 

• 2 9 1113 9 

.9867005 


. 3(36 1 228 

.3013128 

# 301 2556 

.9873925 


.3265306 

.3218700 

#3210102 

.9880398 


.3869388 

.38 t 6809 

• 30 1 5706 

,9887307 


.3673869 

.3618268 

#36)7619 

.9090656 


.3877551 

.3820273 

#3019600 

.9902009 


.8Q8 1633 

.80.22085 

•0021763 

.9910689 


.8285718 

.022079 1 

• 8220090 

.9919380 

P3 f 

.8889796 

.0027318 

.0026609 

.9928523 

Sg 

.8693077 

.0630037 

#06293 l 0 

.9938117 

p0 p 

.8097959 

,0832955 

.0832231 

.9908176 

o a ' 

• 5 10200 | 

.5036085 

• 5035356 

.9958696 


. 5306 1 22 

.5239833 

.5238708 

,9969677 

gg 

.5510200 

.5883010 

# 50022R3 

.9981133 


.57 18786 

.5606826 

#5606100 

.9993062 


.59 10367 

.5850890 

•5850176 

1 . 0005066 

O 

.6127809 

.605521 2 

• 6050509 

1 .001 8351 

p <5 
o Ss 

.6326531 

.6259003 

» 6259 t 1 3 

! .0031722 

§ o 

. 653H6 1 7 

.6068670 

#6063990 

] .0005577 


. 6738698 

. 6 6 69 8 2 6 

#6669170 

1 . 0059928 

MS 

.6938775 

.6875290 

• 6870650 

1 .0070780 


.71 82057 

.7081082 

* 7080030 

1 .0090 1 31 


.7386939 

.7287122 

.7286507 

1 .0-105990 


.755)02" 

*7093531 

• 7092988 

1 .01 22357 


•775510? 

,7700279 

:» 7699770 

1 .0139238 


.7959 100 

.7907377 

• 7906905 

1 .0156603 


• 8 1 63265 

. 8 1 1 8836 

• 81 100P3 

1 .01 70570 


.8367307 

.0322666 

*0322275 

1 .019303 1 


. 857102P 

,8530878 

• 853053 1 

1.021 2020 


.877551 0 

.0739883 

• 8 7 ? 9 1 8 1 

1 .023 1557 


. 8979C.9? 

.0908893 

• 8908230 

1.0251606 


.9 183673 

.9157918 

#9157712 

1 .0272280 


.9387755 

.9367770 

#93676 10 

1 .0293069 


.9591837 

.9578059 

#9577955 

1 .03 15222 


.979591P 

.9708799 

: V; *9788707 

1 .0337507 


1 • OO n On"'-> 

1 • OOOOOOO 

1 • 0000000 

1 .0300850 


NUMBER OF I TfPAT I 0N= 1 



t 


98 


f ■ .10000 PH! ■ .toooo+ni 


z 

VELOCITY 

OLDVEL. 

VStaR 

ooooooo 

•ooooooo 

• ooooooo 

.8474300 

0204Q82 

.0172976 

•0173427 

.8477306 

04it8 1 63 

.03960} 2 

•0346914 

.8481824 

06 1 2245 

.0519172 

.0520521 

.8489375 

Ok 1 6327 

.069251 8 

•0694310 

.0499982 

1 020409 

.0866 1 1 0 

.086834 1 

.8513674 

1 > 2 ^ 0 n 

. 1090014 

. 1 042677 

.8530482 

1428571 

.1214293 

• 1 217380 

.8550441 

1632653 

. ! 38901 2 

• 1 392514 

.8573595 

1336735 

.1564236 

. 1568 1 42 

. 8599986 

7.0408 1 6 

. 1740032 

• 174433 1 

.8629666 

2244898 

. 19 16467 

.1921147 

.8662692 

2848980 

.209361 1 

* 2098656 

.8699126 

2653061 

. 227 1 534 

*2276929 

.8739036 

2fe57 1 43 

.2450306 

.2456036 

.8782492 

306 1 228 

.2630003 

• 2636097 

.8929572 

3265306 

. 28 1 0697 

.28 1 7038 

• 8880364 

3469388 

.2992467 

• 2999083 

.8934958 

3673469 

. 3 1 75389 

.3182260 

.8993453 

3377551 

.3359546 

.3366649 

.9055953 

3081633 

.3545020 

.3552330 

.9122569 

9 2 8 5 7 ! 4 

.3731896 

.3739388 

.9193424 

9489796 

.3920262 

.3927910 

.9268645 

9693877 

.41 10208 

.411 7985 

.9348370 

4397959 

.4301 828 

.4309706 

.9432753 

5 10209 1 

.4495218 

.4503168 

.9521947 

5306 1 72 

.4690479 

.4698470 

.9616116 

5 £ 1 0204 

.4887713 

. 48957 1 5 

.9715494 

57 19286 

.5087028 

.5095008 

.98201 19 

5918367 

.5288534 

.5296462 

.9930344 

6 ! 22999 

.5492348 

.5500189 

1 .0044338 

632653 ! 

.5698589 

• 57063 1 0 

1.0168331 

6 530612 

.5907382 

.5914949 

1.0296569 

6739694 

.6 1 1 8857 

•6126236 

1.0431312 

6938775 

.6333150 

•6340307 

1 .0572840 

7 192857 

.6550401 

.6557301 

1 .0721458 

7396939 

.6770760 

.6777368 

1.0877487 

7551020 

.6994380 

.7000663 

1 . 1041259 

7755 102 

.7221 424 

.7227347 

1.1213144 

79591 84 

.7452060 

.7457589 

1 • 1 393536 

3163265 

.7686466 

•7691569 

1 .1582851 

8367397 

.7924829 

.7929474 

1.1781534 

8571928 

.8167345 

•817150Q 

1 . 1990058 

8 7 7 5 5 10 

.8414219 

•8417855 

1 .2208937 

8 9 7 9592 

.8665669 

.8668758 

1.2438729 

9 183673 

.8921922 

.8924437 

1 .2680020 

9387755 

.9183221 

•9185137 

1 .2933443 

9591837 

.9449818 

•9451115 

1.3199678 

97959 1 8 

.9721 983 

•9722640 

1 .3479459 

OOOOOOO 

1 • ooooooo 

1 •ooooooo 

1 .3622829 


number OF IT f RATION® 3 
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Figure B-4 










CHAPTER II 


INTRODUCTION 


The experimental program mentioned in Chapter I is continued in 
Chapter II with slight modifications of the experimental system. The 
manifold system used in the earlier phase of the program had only four 
three-way stopcocks whereas the new manifold system used seven two-way 
stopcocks. This improvement facilitated in collecting the pressure drop 
data across the channel length at seven locations without disturbing the 
experimental setup. Several precautions were taken especially to eliminate 
airbubbles and to minimize the turbulence in the tank to obtain better 
and more accurate readings. Also two new screens (Dutch twill 200 x 200 
and 325 x 325) were used during this phase of the project to investigate 
the pressure drop across the channel length. Three channel depths, 
namely 2.2, 2.85 and 3.5 cm, were used with the flow rates varying from 
345 cc/sec to 65 cc/sec. 

The second chapter of this portion of the report contains the 
theoretical consideration flow patterns in a rectangular channel followed 
by the summary of the experimental results . 
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SECTION 1 

THEORETICAL CONSIDERATIONS 


The fundamental physical laws of conservation of mass, momentum 
and energy, when applied to a continua, in most cases yield sufficient 
equations for evaluation of the flow parameters. These flow parameters 
are velocity, pressure and forces developed at a given region inside the 
continua. For our analysis we chose a control volume fixed in space 
and evaluated the flow parameters in an averaged integrated fashion. 

The choice of the control volume is arbitrary. However, in order to 
extract maximum information it is essential for the boundaries of the 
surface of the control volume to pass through regions where information 
is known and also where it is required. In most problems, several control 
volumes may be necessary for the formulation of a determinate set of 
equations. For clarity in presentation, we express the vector quantities, 
velocity, momentum and forces in their component form using the 
rectangular cartesian system (an appropriate system for our flow 
geometry). Hence, all the equations are scalar in nature. 

1-1 Equations of Motion 

The channel is shown in Fig. (1) and the dimensions are also 
marked on it. The control volume is shown in Figure (2). Figure (2) 
shows the mass flows, momenta and forces acting on the control volume 
as follows : 

1-2 Conservation of Mass 

Let u and v denote the x and y components of the velocity 
vector and p the density of the fluid flowing through the control 
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volume. The control volume is a rectangular, parallelopiped with the 
dimensions of Ax, D and H in the x, y and z directions of the rectangular 
cartesian system. The voids ratio e is introduced to calculate the 
actual area available for flow due to the pressure of the screen. 

From Sketch (a) : 

Mass flow in = Mass flow out 
vpLAxe + upD L - upDL + ~ Ax p D L 
or 


D du 

v = — ~r~ . 

e dx 


( 1 ) 


In arriving at this result, we make the assumption that the velocity 
component u is a function of x only. This is approximately true because 
below the screen the velocity is predominately in the x direction. 

1-3 Conservation of Momentum 

We apply the law of conservation of momentum in component form 
in the x direction. In other words, we have to calculate the net 
efflux of x-momentum and equate it to the sum of the x-forces acting on 
the control volume. Since the forces in general can act in either direction, 
we adopt the sign convention that forces acting to the right are positive, 
whereas those acting to the left are negativ^. 

The forces on the control volume are due to the pressure, p, of 
the fluid and also due to the shear stress, t g , acting on the control 
surface. Evaluation of the shear stress is not obvious; therefore, we have 
to adopt a suitable representation to do so. A conventional method is 
to introduce a friction factor, f, and to evaluate the losses in energy 
due to friction forces, as a fraction of the unit kinetic energy, u 2 /2g. 
Several representations are therefore possible for relating the quantities 
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T i f and u. One such formula for calculating the wall shear stress, 
r , is to take the corresponding force acting opposite the direction 
of the flow, as t q LD, thus the momentum equation becomes: 

EF x = (momentum out - momentum in) x-direction 
pDL - (p + Ax) DL - T q L D = (u + — Ax) 2 pDL 

( 2 ) 

- [u 2 pDL + v u pLAxe] 

since, T q = ~ P ~ U -- . (3) 

By substituting equation (3) in equation (2) and neglecting terms 
containing Ax 2 and the simplified momentum equation we obtain: 

s + S + f §“ 2 - 0 < 4 > 

The purpose of the analysis is to solve for the flow parameters 
u, v and p. Equations (1) and (4) are insufficient since we have three 
unknowns. Therefore, we must obtain one more equation containing any 
or all of the terms (u, v, p) . 

There is no unique way to generate one more independent equation 
containing u, v and p. We will discuss two different methods, outlined 
as follows. 

Method 1 - A successful method of relating p and v is to involve an 
experimentally determined constant, K (or parameter). Thus, the approach 
becomes quasi- analytical. Let us assume that 

K 0 V = p o - p (5) 

The quantity p Q is a reference pressure introduced for ease in rendering 
the equations dimensionless at a later stage. In differentiating equation 
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(5) with respect to the independent variable x, we obtain 


_ = K dv 

dx o dx * 


( 6 ) 


The term ^ can be evaluated from equation (1) by simply obtaining 
the derivative with respect to x; thus, 


dv 

dx 


D dfu 
e dx 2 


(7) 


Substituting equation (7) in equation (5), we get 


dp K I) d 2 u 
dx “ o e ,jx 2 


( 8 ) 


Equation (8) facilitates elimination of the pressure gradient 
term, dp/dx, in equation (4) and thus yields a single ordinary non- 
linear differential equation for u. This equation is written as 

D d 2 u „ du of 2 n 

K r— - PuT--^u z = 0. 


O G 


dx 


dx 2D 


(9) 


If it is possible to solve this equation by using equation (1) and 
(5), we can determine v and p at any x location inside the control 
volume, as the solution for u from equation (9) is already known. 

Alternatively, a totally different method is proposed here 
which obviates the necessity of the introduction of any new constants 
(e.g., K q as in the method just discussed). 

Method 2 

Consider the control volume shown in figure (2). Applying the 
momentum equation in the Y direction, we obtain 


2 f = (momentum out - momentum in) Y-direction 


or 


P q LAx - pLAx + pgLAxD = v 2 pLAx. 


Differentiating with respect to x we obtain 


iR „ 

dx 


2pe v- 


Substituting equation (1) and (7) for v and ^ respectively in (10) yields 

_ iR , 2p 2i <!fa (11) 

dx e dx dx 2 

Now equation (11) for the pressure gradient can be employed in the x 
direction momentum equation to obtain an equation for u, namely. 


- D z d 2 u du du pf 2 _ n 

2p 2 i — p u i — rtr u - o. 

s dx z dx dx 2D 


It is obvious that these two methods give slightly different 
governing equations for u in terms of the geometry (L, D, H, e) and the 
independent variable x. A very useful conclusion can be drawn at this 
early stage if we compare equations (9) and (12) . These equations become 
identical if K q in equation (9) is chosen as being equal to 2pD ~ . This 
suggests that K q is not an absolute constant but dependent on the velocity 
gradient du/d x. However, it will be a constant if the variation of u in 
the direction of x is linear (not a constant). As it will be seen later, 
the flow rate through the channel has a significant effect on K q . 

2-1 Solution for the Equations of Motion 

In the previous section we discussed the formulation of the 
equations of motion. Two different representations for the variation 
of the u component of velocity were obtained. We consider the detailed 
solution for the equation (9) which is rewritten as. 


D d 2 u du pf 2 

K — — j - pu T- - — u - 0 

° e dx dx 2D 
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To keep the solution general, we have to render equation (9) and other 
related equations dimensionless. For achieving this we use the following 
scheme: 


Q * 
u = — ~ r u 
pDL 


v 


0 * 

— v 


pLHe 


x = HZ 


P 0~ P 

K 


pLHe 


Ap* 


ft ^ 

Here Q is the flow rate through the control volume and u , v and Ap 
are the dimensionless velocity components and pressure differences, 
respectively. The quantity Z represents the dimensionless x co-ordinate. 
All other quantities have been defined earlier. 

On substituting these variables into equation (9), we obtain 


d 2 u* j. * du* , „ *2 

5 cpu cf) F u = 0 

dz 2 dz 


(13) 


The terms $ and F are termed the flow and friction analogs respectively, 
and are defined as: 


<t> = 


JL2_ e 

K 0 D L 


F = 


_f_ 

2D 
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The nonlinear equation (13) is solved by means of the Crank- 
Nicholson method on a digital computer. (A detailed listing of the 


program is included in Appendix A.) 

The flow domain is subdivided into equal intervals of length 
Az between two subsequent stations i and i + 1 or i and i-1. The 
solution for the differential equation is thus a set of finite values 
u^ (i 1, N) and u^-»-u as Az->0. 


The derivatives 


d z u* du* 

r and -3 — of (13) are replaced by the following 

dz dz 


equations : 


,2 * 
d u 

dz 2 


du 

dz 


u i+l ~ 2u i + u i-l 
Az 2 


u i+l " u i-l 


(14) 


2Az 


The resulting algebraic equation, derived from equation (13) by substituting 
equation (14) into it, is linearized by an inter lue scheme as 


A 1 u l-l + *2 U 1 + A 3 u i+l = D „ 


(15) 


where 


A, = 


Az' 


r V-l 

<f> 

2Az 


where 


^ - _ zL . _ p $ u v_1 

2 Az 2 i 


A_ = 


Az' 


$u 


V-l 


2Az 
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V = no of iteration 
V-l = previous iteration. 

The boundary conditions are u^ = 0 at z = 1 and u^ = 1.0 at z = 1.0. 

The calculation is initiated by arbitrarily choosing values for all 
the u^ (v = 1) and then calculating A^, and A^. Since the boundary 

condition at u^ = 1 is known, aback substitution is necessary for evaluating 

V V 

the u^_^ from u^ based on the previous iterant. The procedure is repeated 

_ 4 

until a residue (arbitrarily fixed) of 1 x 10 or less is attained for 
the corresponding improvement in the value for u^ for each progressing 
V. Here we used Az = 10 -2 . 
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3. EXPERIMENTAL RESULTS 

Experimental data with the modified manifold system with seven 
stopcocks were collected for 200 x 200 and 325 x 325 dutch twill screens 
using tap water. The volumetric flow rates were varied from 345 cc/sec 
to 65cc/sec. Three different depths of the channel, namely 3.5, 2.85 
and 2.5 cm were used to collect the data for various flow rates. The 
pressure drops across the screen along the channel length at fourteen 
locations were measured by changing the manifold system only once. 

Figures 3 to 8 present the data on pressure drop across the screen 
along the channel length as a function of dimensionless length z for 
various flow rates and channel depths of 3.5, 2.85 and 2.2 cm. 

Examination of figures 3 to R show the following characteristics of the 
pressure drop across the screen length. 

1. The pressure along the channel length across the screen is 
dependent on the volumetric flow rates. At smaller flow rates the 
pressure drop between the two extreme ends of the channel length is 
small whereas at large flow rates the pressure drop between the two 
extreme ends of the channel length are largest. In otherwords, larger 
the flow rate, larger the pressure drop between the extreme ends of 
the channel length. 

2. Pressure drop across is the Screen along the channel length is 
also a function of the channel depth. At the smallest depth 2.2 cm, 
the pressure drop between the two extreme ends of the channel is large 
compared to the pressure drop observed between the extreme ends of the 
channel at larger channel depths of 2.85 and 3.5. 
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The experimental results obtained in the second phase of the 
program are an improvement over the data collected in the first phase 
of the program and are in general agreement with the conclusions drawn 
in the first phase of the program. 
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APPENDIX 


Computer Program 

Using Crank Nicolson type iteration: 


d 2 u u i+l - 2 u i + "l-l 


dz 2 


Az' 


du* = u i+l ~ u i-l 
dz 


2Az 


The equation becomes 


V „ V . V 

u i+l ~ 2u i + u i-l 

At 2 

Or, 


- $ u 


V V 

v-1 u i+l " u i-l 


2At 


- <1>F u°. uV _1 = 0 
1 l 


, V , . V , . V _■ 

A, u. . + A_ u. + A„ u.j, = D 
1 l-l 2 l 3 i+1 n 


Computational Alogrithem 


* V-1 
1 $U i 

A = — 1 — + 

1 ^ 


A 2 = 


A 3 = 


Az 

1 


- F $u 


V-1 


$ u 


V-1 


Az' 


2Az 


D = 0.0 
n 


V -*■ no. of iteration 
V-1 previous iteration 
Iteration Residue = 1.0 x 10 
In the program Az = 1/100 = 0.01 


i-4 
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Then, 


u 


= 9s*. 

pDL 


_ Qv 

v = 


pLHe 


dAp* _ dv* 
dz dz 


Substituting these variables we get 


V* - - 4p* 

dz 


. *„* id . w u 2* = o 

dz 2 dz 


U) 


* = 


M. 


K D 2 L 


F = 


2D 


From method 2 we obtain. 


d 2 u* du* 
dz dz 


. * du 

$ u * 

dz 


$F u* 2 = 0 


$ = 


eH- 

2D 2 


v* = 


du 

dz 


F = — 
2D 


( 2 ) 


4. (P n ~P) 

Ap _ p l 2 H 2 e = v* 2 = 


du 

dz 


* 2 


If model (2) of the analysis is adopted, we can estimate K q from 
the following equations, 

• K = 20D — = - 2 - — = — 

o M dx H pDL dz LH dz 

For a given q, $ can be evaluated from 

' g|J.‘3 

$ for that $ and some constant f the corresponding 

2D 2 
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dll * 

— , or simply Ap , can be read from the graph and therefore K q can be 
computed theoretically. This can be compared with the value obtained 
from earlier type of calculation. 

Hence to compare Ap from the program (Equation 1) (1) calculate 

x • po ” p qIjHg 

z - — for each tap location. Then Ap*=— (r— — ) , K is experimentally 

n K. u o 

o 

determined. 

Method 2 Procedure 

The equations are 


d 2 u* du* , * du 

$ — 


dz' 


dz 


dz 


- $F u* =0 


F = 


2D 


$ = 


£n 

2D 2 


Here u is independent of Q for a constant ’ f*. 

In turbulent flow, ' f’ remains constant. Therefore for high flow rates 
a channel of given dimensions have a unique flow distribution. However, 

Ap* will be different since, Ap* does involve Q in calculations. 

Ap* = pL 2 H 2 £ 

Q 

Figures 9, 10, 11 and 12 show the variation of u* and Ap* along 
z the non-dimensional variables for two sets of f and cj>. Their effect 
is obvious as only the extreme values are chosen. Also plots for Ap 
at z = 0 and z = 1 into f as the parameter are plotted with respect 
to the flow analog cp. The effect of increasing f is opposite on Ap 
at these two locations. 
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Experimental Verification 

Since only the pressure is measured the most useful graph would 
be that shows the variation of (p q - p) the dimensionless measured 
difference on the manometer taps. 

Procedure 

(1) Calculate K as follows: 

' o 


K = 
o 


( Pq-P> 

ip* 


.pLHSv 

' Q ' 


here (p - p) <- measured pressure drop 

Ap calculated from computer program on read off the graph, 
p, L, H, e, Q are parameters known. 


(2) Hence, for a given screen e and flow rate Q, K q can be found. 


Also 


or 


K V = p - p 
o o 


V - 


P o ~ P 


Therefore V can be calculated. 

Validity of the analysis can be established as follows: 

(a) For small (f) and f, u* V z is a st. line 

s 

du* „ . du , 

u — — = slope is constant, or --- is constant, 
dz dx 


V = — —• = constant from above. Hence for a series of Q’s 
e dx ^ 

small but slightly different each other the test for the constant of v 
proves that the analysis is good for that range. It should start 
deviating after some Q because du*/dz is not linear 0 eq. at <j> = 10.1, 
f = 0.6 as shown by the figures. 

(3) These calculations should reveal the nature of K variation. 

o 
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(4) If model (2) of the analysis is adopted, we can estimate K q from 
the following equations: 

„ „ _ du 2pD Q du* _ 20 du* 

' C o 2pD dx H pDL dz LH dz ' 


For a given Q, 41 can be evaluated from (j) = — r for that <J> and some constant 

2D 


f the corresponding — or simply Ap can be read and therefore K q can 


be computed theoretically and compared with the value obtained 
from earlier types of calculations. 
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R=346 cc/sec 


R=290cc/sec 


R=234 cc/sec 


R= 179 cc/sec 


R= I23cc/sec 


R= 65 cc/sec 



200 X 200 Dutch Twill Screen 
D=2.2cm 



R=346cc/sec 


Q=290cc/sec 

Q=l79cc/sec 


O 


□ Q = l23cc/sec 
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200 X 200 Dutch T\ 
D = 2.85cm 



0 = 346cc/sec 


Q=290cc/sec 


Q=234cc/sec 


Q=65cc/sec 


Z 



6 o 


325 X 325 Dutch Twill Screen 
D = 2. 2 cm 

/ ° 

> Q = 346cc/sec 

;■ / 

/ ° 

3 0=290 cc/sec 

o / 




_/D 

> 0=179 cc/sec 

] 0=123 cc/sec 



0= 65 cc/sec 




Velocity Distribution 
Along Length of Channel 


A => F=O.I , <t>=O.I 
O => F=0.6, 4 )= IO. I 









